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We present molecular dynamics (MD) simulations results for dense fluids of ultrasoft, fully- 
penetrable particles. These are a binary mixture and a polydisperse system of particles interacting 
via the generalized exponential model, which is known to yield cluster crystal phases for the corre- 
sponding monodisperse systems. Because of the dispersity in the particle size, the systems investi- 
gated in this work do not crystallize and form disordered cluster phases. The clustering transition 
appears as a smooth crossover to a regime in which particles are mostly located in clusters, isolated 
particles being infrequent. The analysis of the internal cluster structure reveals microsegregation of 
the big and small particles, with a strong homo-coordination in the binary mixture. Upon further 
lowering the temperature below the clustering transition, the motion of the clusters' centers-of-mass 
slows down dramatically, giving way to a cluster glass transition. In the cluster glass, the diffu- 
sivities remain finite and display an activated temperature dependence, indicating that relaxation 
in the cluster glass occurs via particle hopping in a nearly arrested matrix of clusters. Finally we 
discuss the influence of the microscopic dynamics on the transport properties by comparing the MD 
results with Monte Carlo simulations. 



I. INTRODUCTION 

A common strategy to facilitate the study of the physi- 
cal properties of complex macromolecular aggregates is to 
coarse-grain the intramolecular degrees of freedom [TJ [2] . 
By using standard statistical mechanical tools, it is possi- 
ble to represent each macromolecule as a single point par- 
ticle and to obtain an effective pair potential accounting 
for the free energy of interaction between two such macro- 
molecules. For several macromolecular architectures, in- 
cluding linear chains [3J 0] , rings [51 [3] , stars [7HS] , den- 
drimers [10HI2] or microgels [Ml [14] , the so-derived effec- 
tive potentials are "ultrasoft", i.e., the centers of mass of 
the macromolecules can coincide at a modest energetic 
cost (of order /cbT) without violating excluded volume 
interactions between monomers. 

Ultrasoft particles exhibit a more complex phase be- 
havior than that of hard ones. This originates from the 
interplay between entropy, which governs the structural 
properties of hard-sphere solutions, and energetic contri- 
butions arising from the fully-penetrable character of the 
ultrasoft particles. The topology of the phase diagram 
in the temperature-density plane can be classified in two 
classes: reentrant or monotonic behavior. The behav- 
ior of the Fourier components of the ultrasoft bounded 
potential provides a necessary and sufficient condition 
for observing one or the other class [T5J 116] . Thus, if 
all the Fourier components are positive the crystalliza- 
tion lines are reentrant. A complex cascade of crystalline 
phases is found on increasing the density and these de- 
pend on the specific ultrasoft potential [T7H20] . On the 
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contrary, if the Fourier transform of the potential shows 
negative values the crystallization lines are monotonic in 
the temperature-density plane. The corresponding crys- 
talline phases of this class are non-conventional: the ul- 
trasoft particles form a cluster crystal [21) . This crystal 
consists of clusters of particles located in the nodes of the 
lattice. Another particular feature of this phase is that 
the lattice constant is density independent. A direct con- 
sequence of this property is that the cluster population 
is directly proportional to the density of the fluid [T5] . 

The generalized exponential model (GEM) [21; is a 
well-known example of ultrasoft bounded potential lead- 
ing to the two former scenarios, depending on the spe- 
cific parameters of the model (see Section]!!]). The cluster 
crystal scenario of the GEM has been confirmed in a se- 
ries of computational investigations [TBI I21US5] . Detailed 
investigations of the phase behavior have revealed an ex- 
tremely complex map of cluster crystal structures |26| . 
Some of these works [2"2"H2"5j have focused on the dy- 
namic aspects of cluster crystals, revealing interesting 
properties. The stability of the lattice, which has a non- 
integer average cluster population, is maintained by in- 
cessant hopping of all the particles between the clusters. 
In contrast to the usual observation in glass-forming liq- 
uids [27l [28] , a comparison between Newtonian, Brown- 
ian and Monte Carlo (MC) simulations reveals a signif- 
icant role of the microscopic dynamics on the long-time 
dynamics |25j . In particular the hopping dynamics in 
Brownian and Monte Carlo simulations is characterized 
by short-range jumps, and the long-range, highly direc- 
tional jumps found in Newtonian dynamics are strongly 
suppressed. Recent non-equilibrium simulations of clus- 
ter crystals reveal novel features for their rhcological re- 
sponse |29j . 
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It is worth mentioning that potentials of the cluster- 
crystal class have been derived for some macro- 
molecules [51 [3U] • However, though a certain degree of 
clustering was found in concentrate solutions of such 
macromolecules, cluster crystal formation has not been 
observed yet [5j [31]. The reason is that, because of 
increasing many-body effects at high densities, the ob- 
tained effective potentials were no longer valid at the den- 
sities for which crystallization was predicted. Whether 
one can design specific macromolecules that can form 
cluster crystal phases at high concentrations is still an 
open issue. 

As usual in colloidal systems crystallization may be 
avoided in some situations, leading to amorphous states 
of ultrasoft particles. This allows to investigate the for- 
mation of glassy states for sufficiently high densities or 
low temperatures. To the best of our knowledge, all in- 
vestigations on this issue have been performed in systems 
of ultrasoft particles showing the reentrant scenario for 
crystallization. The counterpart of this phenomenon in 
the amorphous case is a reentrant glass transition. This 
feature was predicted by the Mode Coupling Theory in 
a coarse-grained model of star polymers [33] , though ar- 
rested states could not be investigated because of crystal- 
lization. A natural way of preventing crystallization is to 
introduce dispersity in the particle size. This allowed 
to investigate glassy states of Hertzian spheres [53] . 
confirming the presence of a reentrant glass transition. 
Interestingly, it has been recently shown that Gaussian 
spheres can form glasses at high density even in the ab- 
sence of size dispersity [31] and that these have a strong 
mean- field character |35) . 

In this work we aim to get further insight in the dy- 
namics of ultrasoft particles by investigating the glassy 
behavior for the class of cluster crystal-forming systems. 
We present extensive computer simulations of a binary 
mixture and a polydisperse system of ultrasoft particles 
interacting through the generalized exponential poten- 
tial. We investigate structural and dynamic properties 
around and below the clustering transition. The intro- 
duced size dispersity is sufficient to prevent crystalliza- 
tion and to produce a disordered arrangement of the clus- 
ters' centers-of-mass. We observe the signatures of an 
incoming glass transition, leading to a state that we de- 
note as "cluster glass" , akin to the dynamically arrested 
states observed in colloidal systems with competing inter- 
actions [36 38 . Note, however, that the GEM potential 
is purely repulsive and have no minima, i.e., for the sys- 
tem investigated in this work clustering is found in the 
absence of attractive interactions. A detailed analysis of 
the dynamics reveals a progressive arrest of the clusters' 
centers-of-mass on decreasing temperature, with the re- 
laxation of the particles taking place by hopping between 
the nearly arrested clusters. Finally we provide indica- 
tions that the role of the microscopic dynamics (Newto- 
nian or stochastic) on the long-time dynamics may be 
less important in cluster glasses than in cluster crystals. 

The article is organized as follows. In Section[H]we de- 



scribe the investigated model and give simulation details. 



In Section III we present the simulation results and dis- 



cuss structural and thermo dyna mic properties ( III A.) , as 
well as dynamic properties (III 3). We discuss the depen- 
dence of our results on thermal history and microscopic 
dynamics in subsection |IIip . Conclusions are give in Sec- 
tion rrvi 



II. METHODS 

We investigate the dynamics of ultrasoft fully- 
penetrable particles by means of extensive computer sim- 
ulations of a generalized exponential model of index n 
(GEM-n) [21 . In this model the interactions between 
the two particles are given by the bounded potential 
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For exponents n > 2 the Fourier transform of the poten- 
tial has negative components, and hence the monodis- 
perse system can form cluster crystal phases. In this 
work we focus on two different values of the exponent n 
in Eq. ([TJ: a binary mixture with n = 4 and a polydis- 
perse model with n = 8. 

Binary mixture. — The system is composed of a mix- 
ture of two species {1,2} of particles interacting via a 
GEM-4 potential. The potential is cut and quadratically 
shifted at a distance r c = 2<r a p where a, (3 £ {1,2}. The 
ratio between the particles' diameters is a'22/0'11 = 1.3, 
and the cross-diameter er 12 = (o"n + 022 )/2 = er = 1 is 
set as the unit of length. In this work we focus on the 
case of an equimolar mixture, i.e., with the same number 
of particles for both species 1 and 2. 

Polydisperse model. — The system is composed of N 
polydisperse particles of diameter Uj. Polydispersity is 
introduced by means of a flat distribution of the vari- 
able o~i. The distribution is centered at a = 1 and the 
minimum and maximum values are cr m i n = 0.826 and 
0max = 1-164, respectively. Particles interact through the 
pair potential in Eq. ([I]) with n = 8. The interaction is 
cut off at a distance r c = 1.5<7y, where cry = ((Tj + <Tj)/2, 
CTmin < cry < <7 max , and i,j G {1, ..,N}. In order 
to discriminate between particles of different sizes, we 
introduce three subpopulations of particles labelled by 
a = 1,2, and 3. If we sort the particles by increasing 
value of the diameter <7j, we say that the particle i be- 
longs to the species a if i € [1 + (a — l)A,aA], with 
A = N/3. Thus a increases with increasing average size 
of the particles. 

In both the mixture and the polydisperse system we 
use a common energy scale e a p = 1 and particle mass 
m = 1. The particles are placed in a cubic box with 
periodic boundary conditions. The static and dynamic 
properties of these models were investigated by means 
of Molecular Dynamics (MD) and MC simulations, per- 
formed over a wide range of temperatures T and of den- 
sities p = N/V, with N the number of particles and 
V the volume of the simulation box. Namely we used 
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FIG. 1. Radial distribution functions of the binary mix- 
ture along the isochore p = 4.0 for selected temperatures (see 
legend), (a) ffn(r), (b) gi 2 (r), and (c) 522 (r). 
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FIG. 2. Radial distribution functions of the polydisperse 
model along the isochore p — 5.0 for selected temperatures 
(see legend), (a) gu(r), (b) 013(7-), and (c) g 33 (r). Correla- 
tion functions involving particles of intermediate size (a = 2) 
are not shown. 



N = 4000 in the binary mixture and investigated the 
densities p = 2.0, 3.0, and 4.0. For the polydisperse sys- 
tems we used N = 4394, 4151, 4003 and 2870 for the 
densities p = 2.0, 3.0, 5.0 and 7.0 respectively. 

Production MD runs were performed in the micro- 
canonical ensemble (Newtonian MD). Newton's equa- 
tions of motion were integrated by means of the velocity 
Verlet algorithm [39]. For the polydisperse model the 
time step St ranged from 10~ 3 at high temperature to 
10~ 2 at low temperature. For the binary mixture the 
time step St was 0.02 independent of temperature and 
density. These values of St allowed to keep the degree 
of energy conservation, determined from the ratio be- 
tween the root mean square deviations of the total and 
potential energy, to less than 2% at all the investigated 
state points. Thermalization in the equilibration runs 
was achieved by periodic velocity rescaling in the poly- 
disperse model and by means of the Berendsen thermo- 
stat [40] using a time constant tr = St/Q.l in the binary 
mixture. 

For comparison with the MD results we also performed 
Monte Carlo (MC) simulations in the canonical ensem- 
ble. Propagation of the particles was implemented ac- 
cording to the standard Metropolis algorithm |39J. The 
trial moves performed during the MC simulations in- 
volved random particle displacements, generated over a 
cube of side 0.1. We observed that the acceptance ratio 
varied between a 80% (high T) and a 50% (low T). 

To test the reliability of our results different ther- 
malization criteria were adopted and compared. The 
key quantity to assess equilibration was the root mean 
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evaluated at the end of the simulation. In the case of 
the binary mixture, the duration t cq of the equilibration 
run at any temperature was such that R(t cq ) > 8. The 
duration of corresponding production runs, i pr od, was 
typically four times longer than i cq , so that i?(i pro d) ~ 
2R{t cq ). At the end of the production run, the system 
was cooled to a lower temperature, and the procedure was 
reiterated. To ensure that all the single-particle degrees 
of freedom, i.e., those corresponding to both small and 
large particles, were equilibrated, we also implemented 
an analogous equilibration criterion based on the partial 
RMSD 



Ra(t) = 
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where the sum is restricted to particles of species a. A 
small, systematic dependence on the target value R a (t eq ) 
was observed at the lowest temperatures and will be dis- 
cussed in Sec. MI CI 

In the case of the polydisperse model, equilibration and 
production runs were such that R(t eq ) and R(t pro( i) typ- 
ically exceeded three interparticle diameters. In the runs 
at the lowest investigated temperatures, which covered 
about 10 s time steps, the target RMSD reached values 
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of about one interparticle diameter. Even in these runs 
at very low temperature, however, no drift in the time 
dependence of the potential energy and pressure during 
the production runs could be observed. In addition to the 
same gradual cooling protocol used for the binary mix- 
ture, we also implemented an infinite-quench rate pro- 
tocol, whereby the initial configuration was always pre- 
pared by placing the particles randomly in the simula- 
tion box and then performing an infinite-rate quench to 
the target temperature T, which was subsequently equi- 
librated by monitoring the value of the potential energy 
and pressure as a function of time. We made sure that no 
drift was observed during the production runs. We found 
that the static and dynamic properties of the polydis- 
perse model did not depend appreciably on the quenching 
protocol employed. 

Temperature T, time t, distance r, wave vector q and 
density p are given respectively in units of e/fcs (with ks 
the Boltzmann constant), a(raj e) 1 / 2 , a, er -1 and cr~ 3 . 
Unless otherwise specified, in the following the presented 
results will correspond to the Newtonian MD simulations. 
The comparison between MD and MC results will be dis- 
cussed at the end of Section IIII1 



III. RESULTS AND DISCUSSION 

Unless specified otherwise, in the following we will 
present simulation results for two selected densities: p — 
4.0 for the binary mixture and p — 5.0 for the poly- 
disperse model. As we will show below, the dynamic 
properties of the models investigated here exhibit, at suf- 
ficiently low temperatures, the p/T-scaling found in the 
cluster crystal phase of the monodisperse system [351 122] • 
The scaling becomes effective for p > 3.0 in the binary 
mixture and for p > 5.0 in the polydisperse model (see 
below). The selected isochores are therefore representa- 
tive of the behavior in the high-density scaling regime. 



A. Structure and thermodynamics 

We start our discussion by analyzing the static pair 
correlations. Figures [T] and [2] show the temperature vari- 
ation of the partial radial distribution functions g a /3( r ) 
for the binary mixture and the polydisperse model, re- 
spectively. Both models follow a very similar structural 
evolution along the selected isochores: a prominent peak 
located around r ~ builds up as the temperature de- 
creases, indicating an increased interpenetration of the 
particles. At the same time, the first minimum in <? Q ^(r) 
becomes deeper, suggesting the formation of well-defined 
clusters of typical maximum size tzs 0.7. This value can 
be read off from the positions of the first minima of the 
radial distribution functions. The ability of the particles 
to interpenetrate stems of course from the ultrasoft and 
bounded character of the interaction potential, Eq. ([lj, 
but formation of clusters in purely repulsive models is 
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FIG. 3. Distribution of cluster population numbers rid for 
various temperatures (see legend) in (a) the binary mixture 
and (b) the polydisperse model. 
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FIG. 4. 

(a) for the binary mixture in the (nj^ , 

the polydisperse model in the (n^\n^ 
of the circles are proportional to the probability of finding 
clusters with a given chemical composition. The state points 
are (a) T = 0.35 for the binary mixture and (b) T = 0.45 for 
the polydisperse model. 



a non-trivial collective phenomenon, that arises only at 
sufficiently high density |16j . A close inspection of the 
9ap( r ) m the binary mixture shows that the peaks around 
r rs are significantly higher for like correlations. This 
suggests that in this latter model clusters are mostly pop- 
ulated by particles of the same species — a phenomenon 
that can be described as "chemical segregation". The 
data for the polydisperse model indicate a similar ten- 
dency towards homo-coordination, although the effect is 
significantly weaker than in the binary mixture. 

To characterize these effects more precisely, we per- 
formed a simple cluster analysis. A particle belongs to a 
given cluster if its distance to at least one of the other 
particles of that cluster is smaller than a preselected cut 
off r cu t • When using a fixed value of r cut , it is sometimes 
difficult to unambiguously identify the cluster to which a 
particle belongs, due to the continuous flow of particles 
from cluster to cluster [25] , in particular at high temper- 
ature. In practice, however, we do not observe any major 
artifacts for the systems at hand. In particular, our data 
indicate that merging of neighboring clusters |25j is not 
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FIG. 5. Thermodynamic and cluster properties of the binary 
mixture (black and white symbols) as a function of T : (a) 
total potential energy U(T), (b) specific heat CV(T), and 
(c) fraction P(n c \) of selected cluster populations n c i- The 
vertical dotted line in (b) marks the position of the peak of 
the specific heat. In (a) data for the monodisperse GEM-4 
model at a density p — 4.097 are included for comparison 
(red symbols). 



a serious issue in our models at sufficiently low tempera- 
ture (see Fig. El . In view of the typical widths of the first 
peaks of g a p\r) at low temperature, a reasonable choice 
of the cut-off distance is r cut 0.4 for the binary mix- 
ture and r cut « 0.35 for the polydisperse model. Small 
changes of this parameter (±30%) did not affect qualita- 
tively our analysis. 

Let us first analyze the distribution P(n c \) of cluster 
population numbers n c i. The temperature variation of 
P(n c i) is shown in Fig. [3] for the two models. At high 
temperature the distribution consists of a sharp peak 
around n c \ = 1 (isolated particles), and a featureless 
background of more populated clusters. Thus in the 
high-temperature fluid phase the system is mostly "dis- 
sociated" , though particles may temporarily overlap due 
to collisions. An analysis of the cluster lifetimes (not 
shown) confirms this picture. Moreover, at such high 
temperatures large values of n c i , leading to the observed 
background in P(n c \) 1 may result from the fluid structure 
being spatially more homogeneous, making the definition 
of clusters meaningless. As the temperature is lowered, 
the height of the peak at n c \ = 1 decreases, larger clus- 
ters become more frequent and P(n c \) approaches a well- 
defined ultimate profile. In the polydisperse model at 



FIG. 6. Black and white symbols: as Fig. [5]for the polydis- 
perse model. Red symbols in (a) are data for the monodis- 
perse GEM-8 model at a density p — 5.0. 



p = 5.0, a clear peak is visible around n c \ ~ 8. In the 
binary mixture at p = 4.0 the distribution is clearly bi- 
modal, indicating the formation of two distinct types of 
clusters with average populations of n c \ w 6 and 10, re- 
spectively. We have performed a similar analysis (not 
shown) at the other simulated densities. As expected, 
we find that the "preferred" values of n c \ are density de- 
pendent, namely they tend to increase with increasing 
P- 

A closer inspection of the chemical composition of the 
clusters reveals that the double peak structure of P{n c \) 
in the binary mixture reflects the same chemical seg- 
regation indicated by the radial distribution functions: 
"small" clusters (n c i ~ 6) are mostly formed by small 
particles (species 1), whereas "big" clusters (n c i = 10) 
are mostly formed by big particles (species 2). This effect 
is nicely illustrated in Fig. [4] The population number of 
each cluster is first decomposed as n c ; = +nff , where 

n^p and n? indicate the number of particles of species 
1 and 2 composing the cluster, respectively. Then, the 
two-dimensional map of the distribution of clusters with 
"chemical composition" (tvJ ) is constructed. The 
distributions shown in the figure correspond to thermo- 
dynamic states located in the temperature regime where 
strong clustering is observed. For the polydisperse model, 
the analysis is performed using n^?) and averaging 

over all possible values of n cl . In the binary mixture, the 
tendency towards chemical segregation is rather evident. 
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FIG. 7. Snapshots of the particles' positions, above and 
below the clustering temperature T* of the binary mixture: 
(a) T = 0.75 and (b) T = 0.35. Particles of species 1 and 
2 are depicted as small white spheres and big red spheres, 
respectively. For clarity, only particles contained within a 
vertical slab of thickness 4 are shown. 
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FIG. 8. Same as Fig. [7] but for the polydisperse model: 
(a) T = 2.44 and (b) T = 0.64. Particles of species 1, 2, 
and 3 are depicted as small white spheres, intermediate blue 
spheres, and big red spheres, respectively. For clarity, only 
particles contained within a vertical slab of thickness 4 are 
shown. 



It is worth mentioning that at these low temperatures 
"merged clusters" are extremely rare. Cluster merging 
in the binary mixture would lead to the appearance of 
chemical compositions such as (n^ , rv^ ) = (12,0) or 
(0, 20) (i.e., union of preferred clusters composed of small 
and big particles, respectively), or (6,10) (i.e., union of 
two different preferred clusters). However, the distri- 
bution shown in Fig. [4^i has no contributions around 
,n*j ) = (6,10) and (0,20). Only a few instances 
of (12,0) and (12,1) clusters are visible, the latter aris- 
ing from two (6,0)-clusters merged by one big particle, 
but their fraction is negligible (< 10~ 5 ). Thus, we con- 
clude that the binary mixture of GEM-4 particles self- 



assemble at low temperatures into a binary mixture of 
clusters. The results for the polydisperse model indicate 
an anti-correlation between subpopulations of small and 
large particles, but the effect of chemical segregation is 
much weaker than in the binary mixture (see Fig. |4Jd) . 
Thus, the clusters in the polydisperse model remain in- 
trinsically polydisperse in character. 

In the following, the crossover from the dissociated 
fluid phase at high temperature to the cluster-dominated 
regime at low temperature will be identified, without 
too much rigor, as a "clustering transition" . Though 
we have no evidence of the possible thermodynamic na- 
ture of this phenomenon, inspection of the excess spe- 



7 




FIG. 9. Radial distribution function of the clusters' centers 
of mass in the binary mixture at some selected temperatures 
(see legend): (a) g?^(r), (b) <?f 2 M (r), and (c) gg M (r). 



cific heat CV(T) — j^^f indicates a clear operational 
definition of the transition. The total potential energy 
per particle, U(T), and Cy(T) are shown in Figs. [5] 
and [6] for the binary mixture and polydisperse system, 
respectively. Both systems display a smooth decay in 
U(T) on decreasing T, and a broad peak in CV(T). The 
maximum of this peak, at temperature T* , marks the 
clustering transition. T* takes value 0.52 and 1.3 for 
the binary mixture and polydisperse model, respectively. 
An analogous definition of the clustering transition has 
been used in the study of the microphase separation of 
a model with short-range attraction and long-range re- 
pulsion |41[ [42 . For comparison, we have included in 
Figs. |5j-a and [6 a the respective results for the monodis- 
perse systems (GEM-4 and GEM-8, respectively). These 
systems were prepared in their equilibrium cluster crystal 
phases and subsequently heated up to high temperatures. 
The density was p = 5.0 and p = 4.097 for the GEM-8 
and GEM-4 models, respectively. The latter value of the 
density was chosen so as to match precisely the effec- 
tive packing of the binary mixture within an "effective 
one-component" description [43 . This adjustment was 
necessary to achieve the expected full collapse of the po- 
tential energies of the two GEM-4 models at high temper- 
atures. The abrupt change in the potential energy of the 
monodisperse systems indicates the transition from the 
fluid to the cluster crystal phase, which is a true thermo- 
dynamic transition. In the following we will denote the 
corresponding melting temperature as T m . Interestingly, 
T* of the polydisperse system is very close to the cor- 
responding T m of the monodisperse model, whereas the 
difference is more pronounced for the binary mixture. 

To correlate the variation of the former thermody- 
namic quantities to the formation of clusters in the binary 



FIG. 10. As Fig. [9] for the polydisperse system: (a) 3n M (r), 
(b) 5 % M (r), and (cT S g, M (r). 



mixture and polydisperse system, we show in Figs. [5]- 
c and |6[c the T-dependence of the fraction P(n c i) of 
clusters with selected values of n c ;. The fraction of iso- 
lated particles (n c i = 1) decreases rapidly around T*. 
Concomitantly, the fractions of the preferred clusters in- 
crease, but in a rather smooth fashion. Therefore, al- 
though a clear signature of the clustering transition can 
be found in the thermodynamic properties, one should 
bear in mind that it may well represent a crossover and 
not a sharp, thermodynamic transition. Around T* , clus- 
ters and isolated particles are indeed in continuous and 
dynamic exchange — a sort of "chemical equilibrium" pic- 
ture. Figs. [7] and [8] show typical snapshots of the binary 
mixture and polydisperse system above and below the 
clustering transition T*. The densities (p = 4.0 and 5.0 
for the mixture and polydisperse system, respectively) 
are the same for which structural and thermodynamic 
observables have been presented in previous figures. The 
snapshots confirm visually the general structural features 
discussed above. 

As temperature decreases well below the clustering 
temperature T* , the fraction of isolated particles be- 
comes negligible (< 10%) and the system enters in a 
"cluster phase", i.e., a regime where almost all parti- 
cles form tightly bounded clusters and only rare jumps 
allow particles to be transferred from one cluster to an- 
other. What is the structure of the clusters of such a 
cluster phase? To address this point, we analyze the ra- 
dial distribution functions of the clusters' centers 
of mass (CM). The index a indicates here the "chemi- 
cal composition" of the clusters. Namely, a-clusters are 
defined as those composed by a majority of particles of 
the species a. Isolated particles are excluded from this 
analysis. In Figs. [9] and [10] we show, for temperatures 
below the clustering transition, the functions g ( ^{r) of 
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FIG. 11. Static structure factors of the clusters' centers 
of mass in the binary mixture at some selected temperatures 
(see legend): (a) Sg M (k), (b) S? 2 M (k), and (c) Sg M (k). 



FIG. 12. Static structure factors of the clusters' centers of 
mass in the polydisperse system at some selected tempera- 
tures (see legend): (a) S?^(k), (b) S? 3 M (fc), and (c) Sg M (fc). 



the binary mixture and polydisperse model, respectively. 
We find that the cluster structure changes only slightly 
upon cooling, without any major transformation. In par- 
ticular, the cluster structure of the polydisperse model is 
rather insensitive to temperature variation and clearly 
amorphous. The temperature dependence of the radial 
distribution functions of the binary mixture is somewhat 
stronger than that of the polydisperse model. Moreover 
these functions exhibit a more marked splitting of the sec- 
ond peak, especially for 2-2 correlations. The typical dis- 
tances between neighboring clusters can be read off from 
the location of the first peaks and reflect the different 
"chemical" composition: clusters composed of big parti- 
cles tend to have larger distances to their first neighbors. 
In Figs. 11 and 12 we include the corresponding data for 



the static structure factors of the clusters' CMs, S2]f(k). 
The maxima of the different structure factors take mod- 
erate values, S™(k) ~ 2 for correlations between same 
species and S^g l (k) ~ 1 for distinct species. No signa- 
ture of Bragg peaks, which would be present in cluster 
crystals, is found, confirming the amorphous character 
of the arrangement of the clusters' CMs. The significant 
anticorrelation, S^p l (k) ~ —0.5, of distinct species in 
the low-A: regime suggests a certain segregation between 
clusters of big and small particles. 

In summary, from the former analysis we conclude that 
the structure of the clusters' CMs below the clustering 
transition is essentially amorphous, at least down to the 
lowest simulated temperature, but differences are visible 
between the two investigated models (binary mixture and 
polydisperse system). We will study the dynamic char- 
acter of these amorphous cluster phases (whether they 
consist of fluids or glasses of clusters) in the next subsec- 



tion. 



B. Dynamics 

We now turn our attention to the dynamics of the two 
investigated models. In doing so, we will better char- 
acterize the nature of the cluster phases identified in the 
previous section, and we will highlight the differences and 
similarities with respect to the dynamics of ultrasoft par- 
ticles in cluster crystals [35] . 

Let us start with the analysis of the temperature de- 
pendences of the diffusion coefficients. These have been 
extracted from the long time limit of (R^(t))/6t, where 
the partial mean square displacements are defined as: 



<*£(*)> = < at £i*w 



The index i runs over all particles of species a. In Fig. [13] 
we show the species-dependent diffusion coefficients D a 
of the binary mixture (p = 4.0, panel (a)) and the total 
diffusion coefficient D of the polydisperse model [p = 5.0, 
panel (b)). The very high temperature regime is omitted 
for clarity and will be shown in Fig. |23| An Arrhenius 
representation is used to highlight the development of 
slow dynamics upon decreasing the temperature. A first 
portion of the data, covering temperatures for which the 
diffusion coefficients are between s» 1 and « 10~ 2 , can 
be reasonably well described by a Vogel-Fulcher-Tamman 
(VFT) law 



D a (T)~exp -A^/(T-4 a) ) 
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FIG. 14. Diffusivities of the GEM-8 model vs. p/T. Filled 
symbols are data for the monodisperse system (25]. Empty 
symbols are data for the polydisperse system of this work. 
The thick dashed lines indicate Arrhenius-like behavior in the 
cluster crystal and cluster glass of the monodisperse and poly- 
disperse systems respectively. 




FIG. 13. Arrhenius plot of the diffusion coefficients (sym- 
bols), (a): D\ and D2 for the binary mixture, (b): total dif- 
fusion coefficient D for the polydisperse model. The vertical 
dotted lines indicate the location of the clustering transition 
(T*) and cluster glass transition (T 9 ). Dashed lines are fits 
to an VFT law (for T > T g ) and to an Arrhenius law (for 
T < T s , activation energies are indicated). 



The fitted values of the strength parameters depend 
on both species and models under consideration. We 
observe that the clustering temperature T* lies within 
this first portion of data. Thus no dynamic signature 
of the thermodynamically defined clustering transition 
can be evidenced from this representation of the data 
(see also Fig. 23 1. At lower temperatures the diffusion 
coefficients undergo a crossover to a milder temperature 
dependence, which can be well described by an Arrhenius 
law 



£> a (T)~exp -E%>/T 



The activation energies are species-dependent and tend 
to be higher the bigger the particles. We interpret this 
"fragile-to-strong" crossover (i.e., from VFT to Arrhenius 
temperature dependence) as a signature of a change in 
the microscopic transport mechanisms, from the high-T 
collective dynamics to activated, single-particle hopping 
taking place in a nearly frozen cluster structure. We re- 



mark that this crossover takes place below T* , in a regime 
where clustering is nearly complete and the fraction of 
isolated particles Pi is typically lower than 15-20%. 



Figure 14 shows the diffusivities of the polydisperse 
system versus the variable p/T at the different investi- 
gated densities. At low temperatures, data for p > 5.0 
collapse to a common Arrhenius law, i.e., with an acti- 
vation energy E aci oc p. Low-T data for p = 3.0 do not 
collapse though still are close and parallel to the data for 
p > 5.0. This scaling behavior of the diffusivities in the 
cluster phase of polydisperse GEM-8 particles is analo- 
gous to that observed in the fee cluster crystal phase of 
the monodisperse counterpart. The data of the fee clus- 
ter crystal (taken from Ref. [22 ) are included in Fig. 14 
for comparison. In analogy with the observation for the 
cluster crystals (see discussion in Ref. [22]), the Arrhenius 
temperature dependence suggests that the particles per- 
form hopping dynamics on the sites of an almost frozen 
matrix of clusters. Though the same qualitative p/T- 
scaling is found, quantitative differences are observed. In 
particular the activation free energy in amorphous clus- 
ter phases is lower than in the cluster crystal. This can 
be tentatively assigned to structural disorder leading to 
a larger number of available pathways ( "entropic" contri- 
bution) or to a smoother energy landscape probed by the 
single particles in the amorphous cluster phase, in con- 
trast to the ordered structure of local minima separated 
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FIG. 15. Intermediate scattering functions for the binary mixture evaluated at various temperatures (see legend): (a) 
Fs{k ~ 6, t), (b) Fs(k = 5,t), (c) Fn(k = 6,i) , and (d) F22(k = 5,t). The clustering and cluster glass transitions are 
highlighted with bold dashed and bold continuous lines, respectively. 



by high barriers in the cluster crystal )23j . 

The p/T-scaling of the polydisperse system breaks 
down at sufficiently high temperatures, i.e., above 
the fragile-to-strong crossover. This suggests a more 
complex, cooperative transport mechanism above this 
crossover, in analogy with the mildly supercooled regime 
of glass-forming liquids. We will show below that in- 
deed the systems investigated here exhibit, in this regime, 
characteristic dynamic features of that scenario. Finally, 
it is worth mentioning that the data for the polydis- 
perse system and its monodisperse counterpart show a 
perfect overlap in the high-temperature fluid phase, i.e, 
above the crystallization and clustering transition for 
the monodisperse and polydisperse system respectively. 
Thus, the effects of the size dispersity have no relevance 
for the dynamics at such high temperatures. Consistently 
with the observations for the potential energy and specific 
heat (Figs. [5] and [6]) the diffusivities of the polydisperse 
system display a smooth variation around T* , in contrast 
to those of the monodisperse system around T m . All the 
previous qualitative observations are also found for the 
p/T-scaling (obeyed for p > 3.0) of the diffusivities in 
the binary mixture (not shown). 



To characterize the dynamics in more detail, we study 
the temperature evolution of both incoherent and coher- 



ent intermediate scattering functions (Figs. 15 and 



In the case of the binary mixture, we report separately 
the scattering functions for small and large particles, 
whereas an average over all particles is performed for the 
polydisperse model. The selected wave vectors reflect 
the positions of the respective first peaks in the static 
structure factors. For the binary mixture the values are 
k* = 6.0 (for a = 1) and k* — 5.0 (for a = 2), whereas 
for the polydisperse model k* = 5.8. 

We first discuss the incoherent scattering functions 
F s (k,t). These are calculated as 



F?(k,t) = 



1 



£ aq>[*£ •(»*,(*) -f$(0))], 

.7=1 



(2) 



where the sum is performed over the positions, fj, of the 
N a particles of the species a. At high temperature the 
correlation functions decay rapidly to zero in a simple ex- 
ponential fashion. Around and below T* , marked oscilla- 
tions appear at short times t ~ 1, which can be attributed 
to the vibrations of individual particles within the clus- 
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Polydisperse (a) 




FIG. 16. As fig 15 for the polydisperse model: (a) F s (k 
5.8, t) and (b) F(k = 5.8, i). 



ters. Similar features have been observed in the clus- 
ter crystal phases of similar ultrasoft particles [25] , and 
are known to be associated to single-particle vibrational 
modes [H]. As the temperature is further decreased, 
the damping of these oscillations becomes weaker and 
the amplitude gets larger, consistently with an increased 
stability of the clusters. The oscillations are also visi- 



ble in the polydisperse model (see Fig. 16 1), though they 
are less pronounced due to the average over all particles' 
sizes. 

The appearance of oscillations in F s (k,t) is accompa- 
nied by evident signatures of glassy dynamics. A plateau 
in F s (k, t) develops at intermediate times and apparently 
grows in a continuous fashion from zero to finite val- 
ues. Such an increase of the plateau height indicates 
a dynamical slowing down of the continuous type. In- 
terestingly, this feature resembles the type-A transition 
scenario predicted by the Mode Coupling Theory for cer- 
tain classes of glassy systems [451 US] ■ It should be noted 
that the emergence of the plateau occurs clearly above 
T*, i.e., the onset of slow dynamics already takes place 
prior to the thermodynamically defined clustering tran- 
sition, both for the binary mixture and for the polydis- 
perse system. This is consistent with the smoothness of 
the transition (see Figs. [5] and [6|. Concomitant with the 
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FIG. 17. Structural relaxation times (a) n and T2 for the 
binary mixture and (b) r for the polydisperse model as a 
function of 1/T. The dashed lines represents VFT fits. 



observed decrease of the diffusivities observed in Fig. [13] 
the plateau extends over longer time scales and is fol- 
lowed by a slow decay with increasing relaxation time as 
temperature decreases. 

The coherent scattering functions are calculated as 



(p a (k,t)p p {-k,0)) 



(3) 



where p a {k,t) — Ylj=i ex P ' f$ Inspection of 
the coherent functions F(k* ,t) of the binary mixture 
(Fig. 15 >d) and the polydisperse model (Fig. 16 o) re- 



veals similar features of glassy dynamics. Since coher- 
ent functions probe collective correlations, the data of 
Figs.[T5"]>d and reflect a progressive dynamic arrest of 
the cluster structure. At very low temperature, F(k*,t) 
does not relax to zero within our observation times and 
the cluster structure is effectively frozen. The highest 
temperature at which F(k*,t) does not completely re- 
lax to zero over the available time scale defines the clus- 
ter glass transition and will be denoted as T g . We find 
T g s» 0.48 and T g « 1.02 for the binary mixture and the 
polydisperse model, respectively [4"7] , 

The slowing down of F(k*,t) above T g follows the 
two-step scenario typically observed in supercooled liq- 
uids. This allows one to define, as usual, the corre- 
sponding structural relaxation times r from the condition 
F(k*,r) = 1/e. Figure 17 shows the temperature varia- 
tion of the species-dependent relaxation times t\ and T2 
for the binary mixture and the total relaxation time r for 
the polydisperse model. The data have been tentatively 
fitted using the VFT function, yielding values of Tq in 
the range 0.42-0.45 for the binary mixture and 0.65-0.7 
for the polydisperse model. These temperatures can be 
taken as lower bounds for the cluster glass transition, cor- 
responding to infinitely slow cooling rates. Furthermore, 
the separation between T* and To is very pronounced in 
the polydisperse model, which indicates that the cluster 
glass formed by the latter is "stronger", in the Angell's 
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FIG. 18. Typical displacements \fi(t) — fi(to)\ of selected small particles in the binary mixture at T = 0.3843 (panels (a) 
and (b)) and T = 0.4983 (panels (c) and (d)). As indicated in the legend, the color of the line indicates the population of the 
cluster to which the particle belongs at time t. Portions of the trajectory during which the particle is isolated are highlighted 
with a thick red line. 



terminology [48] , than the one formed by the binary mix- 
ture. All this clearly shows that the clustering and clus- 
ter glass transitions are distinct and different in nature: 
the former is an equilibrium crossover, the latter an er- 
godicity breaking transition that would occur at different 
temperatures depending on observation times. 

Let us summarize the scenario emerging from our simu- 
lations. Consistent with the structural features discussed 



in subsection III \, the clustering transition is associated 
to a crossover from a fluid of mostly dissociated parti- 
cles to a fluid of clusters. The dispersity of the clusters 
evidenced in |III^ V frustrates crystallization of the clus- 
ters CMs and allows quenching the fluid to lower tem- 
peratures. In the temperature range T g < T < T* the 
underlying transport processes are rather complex: by vi- 
sual inspection of animated particles' trajectories, we ob- 
served diffusive motion of clusters as a whole, branching 
of clusters as well as random walks of isolated particles. 
This results in non-Arrhenius temperature dependence 
of the transport coefficients and a non-trivial behavior 
of the intermediate scattering functions. Below T g the 
structure of the clusters' centers-of-mass becomes prac- 
tically arrested on the time scale of our simulation (see 
below) and the system enters in the cluster glass regime. 
Having noted this, the incoherent scattering functions 
still decay to zero in the time scale of the simulation at 



ergodic in the cluster glass and is driven by hopping 
between the effectively arrested clusters. This hopping 
dynamics is manifested by the Arrhenius-like behavior 



observed in the diffusivities (see Fig. 13). Therefore, as 



in cluster crystals [22], single-particle and collective dy- 
namics in the cluster glass are strongly decoupled. Rep- 
resentative particles' displacements illustrating the dif- 
ference in transport below and above T g are depicted 
in Fig. [I8^i-b and Fig. 18 >d, respectively. Panels (a) 



some temperatures T < T g (see Figs. 15i-b and 16 a,). In 



other words, single-particle dynamics is still effectively 



and (b) show typical single-particle hopping processes at 
low T, involving jumps over a few neighbor distances. 
Note that, during the jumps, the motion of the particle 
is nearly ballistic, as in cluster crystals [25] . and may in- 
volve passing through different clusters. Panels (c) and 
(d) illustrate the more complex motions occurring above 
T g , which involve both slow diffusion of particles resid- 
ing within clusters (panel (c)) and diffusion of isolated 
particles (panel (d)). 



To conclude this section, we provide now explicit evi- 
dence that the cluster glass transition corresponds to the 
loss of ergodicity in the degrees of freedom associated to 
the clusters' centers of mass. At each time t we identify 
the positions Rj(t) of the clusters' centers of mass, where 
1 < 3 < N c i (t) and N c i (t) is the total number of clusters 
in the system at time t. We then calculate the coherent 
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FIG. 19. Symbols: coherent scattering functions F^^ l (k*,t) evaluated for the clusters' centers of mass in the binary mixture. 
The respective data for the coherent scattering functions F a g(k,t), calculated on a particle basis, are included as solid lines. 
Panels (a) and (b) show data at different temperatures (see legend) and fixed wave vector. The selected wave vectors are k* = 6 
for 1-1 correlations (a) and k* — 5 for 2-2 correlations (b). Panels (c) and (d) show data at fixed temperature T = T g = 0.48 
and different wave vectors (see legend) for 1-1 correlations (c) and 2-2 correlations (d). 



scattering functions of the clusters' CMs as 



F2f{k,t) 



<pS m (M)^ m (- 



•M)>' 



(4) 



where yO^ M (fc, t) = Ylj^i ex P ' Rf(t)] an< ^ the sum is 
done over the N°i(t) clusters of the species a. The chemi- 
cal species a of a given cluster is again defined as the ma- 
jority species in the cluster (see above). Single clusters 
(i.e., isolated particles) are excluded from this analysis. It 
is worth mentioning that the averages in Eq. @ are per- 
formed in a grand canonical ensemble, since the number 
of clusters in the sample fluctuates in time. A compari- 
son between F a p(k,t) and F^} /[ (k,t) of the binary mix- 



ture is shown in Fig. 19 |a-b), for wave vectors at the first 



peak in the corresponding static structure factors, and 
at different temperatures. Both data sets follow closely 
each other. When looking at the former scattering func- 
tions at some fixed, low temperature, but changing the 
wave- vector (Fig. [l9|c-d)) the agreement remains overall 
good, especially between big particles and big clusters 
(2-2 correlations). With this, we conclude that the col- 



lective slowing down of the fluid is indeed driven by the 
progressive arrest of the clusters' structure. 



C. Dependence on thermal history and microscopic 
dynamics 

To corroborate our interpretation of T g as a glass 
transition occurring at the clusters' level, we now study 
how the clusters' structure and dynamics depend on the 
quench history obtained from three different thermal his- 
tories for the binary mixture. The difference between the 
data sets considered in this section lies in the time t cq al- 
lowed for equilibration: the equilibration runs are such 
that the partial RMSD R a (t cq ), evaluated at the end of 
the run, is always larger than one, two, and three, respec- 
tively. Note that this equilibration criterion is slightly 
different from the one employed in the rest of the work 
(indicated in Fig. 20 for comparison), which is based on 
the total RMSD of particles. Production runs are typi- 
cally 4 times longer than the corresponding equilibration 
runs [35] . In Fig. 20 we analyze the influence of the 
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FIG. 20. T-dependence of the diffusion coefficients D\ and 
D2 in the binary mixture for different equilibration criteria. 
The different data sets correspond to different values of the 
total RMSD (open symbols) and partial RMSD (full sym- 
bols) targeted during the equilibration run (see legend). In- 
set: equilibration time t Bq as a function of T for the different 
equilibration criteria. 



different quench rates on the T-dependence of the diffu- 
sion coefficients. We remark that, since the equilibration 
times are adjusted at each thermodynamic state to ob- 
tain the target RMSD, a measure of the average quench 
rate is not appropriate. In fact, as shown in the insets 
of Fig. [20j the equilibration times increase exponentially 
with decreasing T. Within statistics the three data sets 



of Fig. 20 show a perfect overlap for T > T g . However, 



small but discernible differences are visible below T g (er- 
ror bars are of the order of the symbol size). Thus, longer 
equilibration runs lead to slightly lower diffusion coeffi- 
cients. A similar and more pronounced effect is visible 
in the relaxation of the coherent intermediate scattering 
functions, which are shown in Fig. [21] at T = 0.48 = T g 
for two different quench rates. Not only r increases with 
increasing R a (t eq ), but also the extent to which F a (k,t) 
has relaxed during the simulation increases with increas- 
ing R a (t eq ), hence with the total length of the production 
run. This, in turn, shifts the cluster glass transition to 
lower temperatures. 

Whereas at the superficial level these observations re- 
semble the typical aging process in structural glasses, 
the microscopic origin of the quench rate dependence of 
the diffusion coefficients is not trivial. In fact, at least 
two factors may affect the single particle dynamics be- 
low the cluster glass transition. First, the structure of 
the clusters' CM relaxes faster upon decreasing equili- 
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FIG. 21. Dependence of the coherent intermediate scattering 
functions (a) F\(k = 6,t) and (b) i<2(fe = 6,t) on the equi- 
libration criterion in the binary mixture at T = T g = 0.48. 
Colors are the same as in Fig. |20| 
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FIG. 22. Dependence of the radial distribution functions (a) 
5ii M ( r ) an d (b) <?22 M ( r ) on the equilibration criterion in the 
binary mixture at T — 0.3843 < T g . Colors are the same as 
in Fig. [20] 



bration times. The single particle dynamics, in turn, is 
affected by the relaxation of the cluster structure, since 
some relaxation channels (e.g., diffusion of a cluster as a 
whole), which would be suppressed upon longer anneal- 
ing, remain active and increase the diffusivity. 

One may also speculate that a further contribution to 
the single particles' mobility arises from differences in the 
topology of the cluster glass, which determines in turn 
the available pathways for particles' hopping. To assess 
whether this is a possible scenario, we analyze the radial 
distribution functions gf4 M (i?) of the clusters' CM in the 
binary mixture as a function of the equilibration crite- 



rion (see Fig. 22). We see that the cluster structure of 
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the most slowly annealed sample differs appreciably from 
the others. Though the positions of the peaks are un- 
changed, their amplitude is increased, indicating a more 
pronounced ordering. These somewhat large differences 
can be explained by the fact that, as evidenced in the 
inset of Fig. [20j the average quench rate does not change 
linearly with the value of R a (t cq ). |50j From the compari- 
son of cluster crystals and cluster glasses (see Fig. 14 ) , we 



see that hopping in a more disordered matrix decreases 
the activation barrier for diffusion, and hence increases 
the particles' mobility. The actual dependence of D on 
the quench rate shown in Fig. [20] suggests therefore that 
a similar effect might be at play in this model. More sys- 
tematic investigations should be performed to assess how 
the specific structure of the disordered matrix affects the 
transport properties at the single particle level. 

Finally, we briefly analyze the role of the microscopic 
dynamics on the transport properties. Following |25| we 
compare the diffusion coefficients obtained from New- 
tonian and Monte Carlo dynamics for the polydisperse 
model. The results of this analysis in the polydisperse 
system are shown in Fig. 23 For temperatures above the 



clustering transition the diffusion coefficients obtained 
from the two methods differ qualitatively above T*, in 
that the MC data saturate at high temperature, whereas 
the MD data show a monotonic increase. Such differ- 
ences are expected, since at high T the GEM-n inter- 
actions will be negligible (random walk on a nearly flat 
energy landscape). Therefore in that limit the MC dy- 
namics will become purely random and T-independent. 
On the contrary, in Newtonian dynamics the onset time 
of the crossover from ballistic to diffusive motion in the 
mean square displacement increases systematically with 
increasing temperature, and so does the diffusivity. In- 
terestingly, in this regime, the temperature dependence 
of D is approximately given by a power law T v , with 
v « 2.2. This value of v is higher than the ones ex- 
pected for both ideal gas (y — 0.5) and Brownian motion 
(y = 1), and remains to be understood. 

Only below T g (i.e., in the cluster glass) the data sets 
obtained from the two dynamics can be collapsed rea- 
sonably well onto a unique master curve, by appropri- 
ately rescaling the time units. With this, the micro- 
scopic dynamics does not seem to play a relevant role 
in the transport properties of the cluster glass. This re- 
sult is rather different from the observation in cluster 
crystal phases [25]. In such systems the Newtonian dy- 
namics leads to long persistent jumps over the lattice 
sites. These highly directional motions are suppressed 
in the MC dynamics, and the corresponding MD and 
MC diffusivities in the cluster crystal exhibit a qualita- 
tively different T-dependence, i.e, they cannot be mutu- 
ally rescaled by a T- independent time factor 25j. The 
fact that, instead, scaling of Newtonian and MC dynam- 
ics works reasonably well in the cluster glass suggests that 
the disordered structure of the clusters' centers-of-mass 
partly suppresses the mentioned long jumps in Newto- 
nian dynamics. A more systematic analysis is required 
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FIG. 23. Comparison of the temperature dependence of the 
diffusion coefficients obtained from Newtonian (circles) and 
Monte Carlo (squares) dynamics in the polydisperse system. 
The diffusion coefficients of the MC data set have been mul- 
tiplied by a factor 75. The vertical dotted lines indicate the 
clustering transition T* and glass transition T g . The dashed 
lines indicate power-law and linear behavior, and are included 
for comparison with the simulation data. 



to completely settle this issue. Work in this direction is 
in progress. 

IV. CONCLUSIONS 

We have presented a computational investigation of 
the structural, thermodynamic and dynamic features of 
dense fluids of ultrasoft fully-penetrable particles. The 
investigated systems are a binary mixture and a polydis- 
perse system of particles interacting via the generalized 
exponential model, for which equilibrium cluster crystal 
phases exist in the monodisperse case. Because of the in- 
troduced dispersity, the systems investigated in this work 
form disordered cluster phases. 

We have characterized the structure of the clusters. 
The analysis reveals a microsegregation of the big and 
small particles, and a strong homo-coordination in the 
case of the binary mixture. The analysis of the thermo- 
dynamic observables does not provide yet evidence for 
a thermodynamic transition associated to the clustering 
transition. Instead, the clustering transition appears as 
a smooth crossover to a regime in which most of the par- 
ticles are located in clusters, isolated particles being in- 
frequent. The structure of the clusters' centers-of-mass 
mirrors that of the dissociated fluid: the binary mixture 
effectively self-assemble into a binary mixture of clusters, 
whereas in the polydisperse model the clusters remain 
polydisperse in population and size. 

The dispersity of the clusters drives a progressive dy- 
namical arrest at the level of the clusters' centers-of- 
mass on approaching the "cluster glass transition" , T g . 
The latter is operationally defined by inspection of the 
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coherent scattering functions, which exhibit character- 
istic features of glass-forming liquids above T g and do 
not completely relax to zero below T g . The diffusivi- 
ties exhibit a fragile-to-strong crossover upon decreasing 
temperature. The onset of strong (Arrhenius-like) be- 
havior occurs below the clustering transition and signals 
a change in the transport mechanisms. Relaxation be- 
low the cluster glass transition is driven by particle hop- 
ping between the nearly arrested clusters. The analysis 
of the dependence of dynamic and structural properties 
on the equilibration times confirms the glassy nature of 
the systems at low temperatures. Finally, we have in- 
vestigated the role of the microscopic dynamics in the 
transport properties, by comparing results from Newto- 
nian and Monte Carlo simulations. At low temperature, 
the diffusivities obtained by both methods can be related 
reasonably well by a single scaling factor. This suggests 



that the microscopic dynamics might play a less signifi- 
cant role for cluster glasses than for cluster crystals. In 
view of the recent observation of clustering in fully atom- 
istic models of interpenetrable colloidal particles [3T] , we 
believe disordered cluster phases of ultrasoft particles de- 
serve deeper and more systematic investigations. 
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